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Abstract 

For atoms with large scattering length, the dependence of the 3-body recombination rate on 
the collision energy is determined by the scattering length and the Efimov 3-body parameters and 
can be expressed in terms of universal functions of a single scaling variable. We use published 
results on the 3-body recombination rate for 4 He atoms to constrain the universal functions. We 
then use those universal functions to calculate the 3-body recombination rate for other atoms with 
large scattering length at nonzero temperature. The constraints from the He results are strong 
if the scattering length is near the minimum of the 3-body recombination rate at threshold. We 
apply our results to 133 Cs atoms with a large positive scattering length, and compare them with 
experimental results from the Innsbruck group. 
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INTRODUCTION 



Vitaly Efimov discovered in 1970 that 3-body systems of nonrelativistic particles with 
short-range interactions have remarkable universal properties if the S-wave scattering length 
a is large compared to the range. If a = ±00, there are infinitely many 3-body bound states 
(Efimov states or Efimov trimers) with an accumulation point at the scattering threshold 
and a geometric spectrum [l[ : 

EP = {e- 2n/so ) n ~ n *h 2 K 2 Jm, (1) 

where is the binding wavenumber of the branch of Efimov states labeled by n*. This 
geometric spectrum is a signature of a discrete scaling symmetry with discrete scaling factor 
e 7r/so |2|. i n the case of identical bosons, so ~ 1.00624 and the discrete scaling factor is 
gw/so ~ 22.7. Efimov showed that the discrete scaling symmetry is also relevant for finite a 
[3|. There is an infinite sequence of negative values of a differing by e n ^ s ° and approaching 
—00 for which there is an Efimov trimer at the 3-atom scattering threshold. There is also 
an infinite sequence of positive values of a differing by e 71 ^ and approaching +00 for which 
there is an Efimov trimer at the atom-dimer scattering threshold. We will refer to the 
universal phenomena characterized by a discrete scaling symmetry that can occur in 3-body 
systems with a large scattering length as Efimov physics. Another example of Efimov physics 
that was discovered more recently is an infinite sequence of positive values of a differing by 
e 7r/so anc j approaching +00 for which the 3-body recombination rate into the shallow dimer 
vanishes at the threshold j3, S, @]: a = (e 7r//s °) n a*o, where a*o ~ 0.32 k~ x . 

In some systems, such as 4 He atoms, the universal aspects of Efimov physics are deter- 
mined by two parameters: the scattering length a and the Efimov parameter k*. However, 
the situation is more complicated for the alkali atoms that are used in most cold atom exper- 
iments as they form diatomic molecules with many deep- lying bound states (deep dimers). 
Efimov physics is modified by the existence of the deep dimers, because Efimov trimers can 
decay into an atom and a deep dimer. The deep dimers also provide additional 3-body 
recombination channels. If there are deep dimers, the universal aspects of Efimov physics 
are determined by three parameters: a, k*, and a parameter 77* that determines the widths 
of Efimov trimers 0] ■ 

The first experimental evidence for Efimov physics has recently been presented by the 
Innsbruck group |8|. They carried out experiments with ultracold 133 Cs atoms in the lowest 
hyperfine state, using a magnetic field to control their scattering length. They observed 
a resonant enhancement in the 3-body recombination rate at a ~ —850 a that can be 
attributed to an Efimov trimer near the 3-atom threshold. At the temperature 10 nK, 
the loss rate as a function of a can be fit rather well by the universal formula for zero 
temperature derived in Ref. [7J with a width parameter 77* = 0.06(1). The Innsbruck group 
also observed a local minimum in the 3-body recombination rate near a ~ 210 ao that might 
be attributable to Efimov physics. One complication is that a is not large compared to 
the van der Waals length scale for Cs atoms: (mCg/a 2 ) 1 / 4 fa 200 ao- Universal predictions 
might not be quantitatively accurate for such a small value of a. Another complication is 
that the measurements were carried out at 200 nK. At this temperature, there could be 
large corrections to the universal predictions for zero temperature. 

The 3-body recombination rate can be calculated at nonzero temperature by carrying out 
a thermal average of the 3-body recombination rate as a function of the collision energy. The 
universal predictions for nonzero collision energy have not yet been calculated. There have 
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been several previous efforts to calculate the 3-body recombination rate for atoms with large 
scattering length at nonzero temperature. D'Incao, Suno, and Esry calculated the 3-body 
recombination rate at nonzero temperature for a simple model potential whose range and 
depth were tuned to give a single S-wave bound state and a large scattering length [§]. They 
considered the qualitative effect of temperature on both the resonant enhancement in the case 
a < and on the local minimum in the case a < 0. Jonsell [l(| and Yamashita, Frederico, 
and Tomio [11] considered the effects of temperature on the resonant enhancement in the 
case a < 0. Jonsell used the adiabatic hyperspherical approximation to calculate the 3-body 
elastic scattering rate as a function of the collision energy Yamashita et al. calculated 
the position and width of the Efimov resonance using a numerically exact method (llj . 
The two groups used their results to estimate the 3-body recombination rate at nonzero 
temperature. They were able to describe qualitatively the temperature dependence of the 
data from the Innsbruck experiment. Neither of these calculations is a definitive universal 
prediction for the 3-body recombination rate at nonzero temperature. 

In this paper, we use published results on the 3-body recombination rate for 4 He atoms 
to constrain the universal functions that determine its dependence on the collision energy. 
Those universal functions are then used to calculate the 3-body recombination rate at 
nonzero temperature for other atoms with large positive scattering length. The constraints 
from the 4 He results are strong if the scattering length is near the minimum of the 3-body 
recombination rate at threshold. We apply our results to 133 Cs atoms with a large posi- 
tive scattering length, and compare them with the experimental results from the Innsbruck 
group. 

THREE-BODY RECOMBINATION 

Three-body recombination is a 3-atom collision process in which two of the atoms bind to 
form a diatomic molecule (a dimer). The 3-body recombination rate R is a function of the 
momenta p 1; P2, and P3 of the three incoming atoms. Galilean invariance implies that R 
does not depend on the total momentum P tot = Pi +P2 +P3- It can therefore be expressed 
as a function of a pair of Jacobi momenta P12 and P3,i2- It is convenient to parameterize the 
Jacobi momenta in terms of the collision energy E, four angles describing the orientations 
of the Jacobi momenta, and a hyperangle a 3 determined by the ratio of the magnitudes of 
the Jacobi momenta. We denote the 3-body recombination rate averaged over the angles 
and the hyperangle by K 3 {E): 



If the scattering length a is positive and large compared to the range of the interaction, 
one of the dimers that can be produced by the recombination process is the shallow dimer 
with binding energy 



If there are deep dimers, they can be produced by the recombination process for either sign 
of a. The recombination rate can be decomposed into the contribution from the shallow 
dimer and the sum of the contributions from all the deep dimers: 



K 3 (E) = (#(pi2,P3,12))p 



Pl2,P3,12, Ct3 ■ 



(2) 




(3) 



K 3 (E) = K shallow {E) + K deep (E) . 



(4) 
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The 3-body recombination rate into the shallow dimer can be further decomposed into 
contributions from the channels in which the total orbital angular momentum of the three 
atoms has definite quantum number J: 



^shallow(£) = $^ (J) (£)- (5) 
J=0 

The scaling behavior near the scattering threshold for each of the angular momentum con- 



tributions is known |l2|: K^(E) ~ E Xj , where Aj = 0, 3, 2, and 3 for J = 0, 1, 2, and 3. 
At the scattering threshold E = 0, only the J = term is nonzero. 

In a gas of atoms with number density tia, the rate of decrease in the number density 
due to 3-body recombination defines an experimentally measurable loss rate constant L 3 : 

—n A = -L 3 n\ . (6) 

The event rate constant a for 3-body recombination is defined so that the number of re- 
combination events per time and per volume is otn\. If the number of atoms lost from the 
system per recombination event is ni ost , the rate constant is L 3 = ni ost a. The binding energy 
of the dimer is released through the kinetic energies of the recoiling atom and dimer. If their 
kinetic energies are large enough that the atom and dimer both escape from the system and 
if they don't undergo further collisions before escaping, then rii ost = 3. On the other hand, 
if they both remain in the system and if we regard the dimer as a distinct chemical species, 
then niost = 2. 

In an ensemble of identical bosons with number density and 3-particle momentum 
distribution /( Pl , p 2 , p 3 ), the event rate per volume and per time is 

an A = - J R( Pl2 , p 3 , 12 ) /( Pl , P2 , P3 ) J^J^J^y 3 ■ (7) 

The factor of 1 /6 compensates for overcounting the states of the three identical particles by 
integrating over the entire range of the three momenta. The 3-particle momentum distribu- 
tion is normalized so that the number of 3-particle states per volume-cubed is n\/?>\\ 

1 f fl x d 3 p! d 3 p 2 d 3 p 3 _ n\ 

6 J /(Pl ' P2 ' P3) (27r)3(27r)3(27r)3-T- (8) 

Dividing Eq. © by Eq. (jSJ), we get an expression for a as a ratio of two integrals. 

If the system is in thermal equilibrium at temperature T, the 3-particle momentum distri- 
bution /(pi , p 2 , P 3 ) is the product of three Bose-Einstein distributions. If T is large com- 
pared to the critical temperature for Bose-Einstein condensation, ksT c ^> 3.3 t^n 2 / 3 /m, the 
3-particle momentum distribution can be approximated by a Boltzmann distribution propor- 
tional to exp(—E tot /(kBT)), where E tot is the total energy of the three atoms. The total en- 
ergy is the sum of the center-of-mass energy and the collision energy E: E tot = P t „ t / (6m)+E. 
The Boltzmann factor is therefore the product of a Boltzmann factor for the center-of-mass 
energy and a Boltzmann factor for the collision energy. Since the recombination rate does not 
depend on the total momentum P tot , the event rate constant a(T) reduces to the Boltzmann 
average over the collision energy E: 

, rdEE 2 e- E ^ k ^K 3 (E) 



6f™dEE*e- 
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The weight factor E 2 comes from using hyper spherical variables for the Jacobi momenta. 
The integral in the denominator can be evaluated analytically to give 2(/c#T) 3 . 



UNIVERSAL BEHAVIOR 

In this section and in the subsequent two sections, we consider atoms for which the only 
diatomic molecule is the shallow dimer and there are no deep dimers. The total 3-body 
recombination rate K^{E) is therefore equal to K s h a \\ ow (E) . An example of an atom with a 
large positive scattering length and no deep dimers is 4 He. In this section, we summarize 
the universal information that is known on the recombination rate Ks(E) for this case. 

At the scattering threshold E — 0, only the J = term in the decomposition of 
-^shallow {E) into orbital angular momentum channels in Eq. ([5]) is nonzero. An analytic 
expression for this term has recently been derived [l3|,[l4|: 



K<°\E - 0) - 7687r2 ( 4?r ~ sin 2 [s In (a/q, )] ha 4 
sinh 2 (7rs ) + cos 2 [s ln(a/a* )] m 

The parameter a„, differs from ft" 1 by a multiplicative constant that is known only to a 
couple digits of accuracy: a* f& 0.32 k~ x . The coefficient of ha 4 /m in Eq. (fTDl) is a log- 
periodic function of a/a*o that oscillates between zero and a maximum value C max as a 
function of ln(a). The analytic expression for the maximum value is 1 

7687r 2 (47r - 3y/3) 
sinn {ttso) 

Its numerical value is C max ~ 402.7. The expression in Eq. fflUj) has zeroes when a is 
(e 71 "^ )™ a*o, where n is an integer. Since sinh 2 (7rso) ~ 139 is so large, K^(0) can be 
approximated with an error of less than 1% of C max /ia 4 /m by 

K®\E = 0) « C max sin 2 [s ln(a/a* )] Tia 4 /m. (12) 

This approximate functional form of the rate constant was first deduced in Refs. 0, 0]. 
The coefficient C mSuX and the relation between a*o and «;* was calculated numerically in 
Refs. 0, [3]. The maxima of the coefficient of ha 4 /m in K(°\0) in Eq. ffl2|) occur when a is 
(e 7r / So ) n 4.76 a*o- The maxima of K^(0) in Eq. (1121) occur when a is (e 71 "/ 50 )™ 14.3 g^q- 

If the collision energy E is small compared to the natural energy scale h 2 /(m£ 2 ) set by 
the range £, the 3-body recombination rate X s h a iiow(-E') is a universal function of the collision 
energy E, the scattering length a, and the Efimov parameter a^. The calculation of the 
rate can be reduced to the calculation of universal functions of a scaling variable x defined 
by 

x= {ma 2 E/h 2 ) 1/2 . (13) 

Only the J = term in A^ shallow (i?) depends on a* , and that dependence is log-periodic with 
discrete scaling factor ~ 22.7. The constraints of universality are therefore particularly 
simple for J > 1. These terms are determined by a single universal function for each J: 

K {J \E) = fj{x)ha 4 /m. (14) 



This constant differs from the constant C max in Ref. 



15( by a factor of 6. 
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The constraints of universality are more intricate for J = 0. They can be deduced from 
Efimov's radial law, which follows from unitarity and from the fact that there is only one 
attractive adiabatic hyperspherical potential in the scale-invariant region of the hyperradius 
R: I <C R <C a. Efimov's radial law implies that the J = term in the 3-body recombination 
rate must have the form 0] 



K^\E) 



x' 



S3i{x)s 12 {x)e 



2iso ln(a/a t o) 



m 



(15) 



1 - Sll (x)e 2isoln ( a / a *°) 

where k is a constant and the functions Sy(sc) are entries of a symmetric 3x3 unitary matrix 
The values of s u (x) and Si 2 (x) at the threshold are known analytically up to phases [T" 



sn(0) 
si 2 (0) 



g-47rso e 2 *^ 2 



(16a) 
(16b) 



The entries s^ix) and ssi(x) go to zero as x 2 as the threshold is approached. Their ratio 
at the threshold is also known analytically (lij : 

(17) 



s 3 x{x)/s 32 {x) — > ^coth(ns )e 2i{Sl - S2) as x -> . 



This limiting behavior together with the analytic values in Eqs. f fT6l) are sufficient to deduce 
the analytic result in Eq. ffTUj) up to a multiplicative normalization constant provided that 
St = ±tt/2. 



ANALYSIS OF 4 HE ATOMS 

4 He atoms are a classic example of atoms with a large positive scattering length. There 
are several modern potentials that are believed to describe the interactions of 4 He atoms 
with low energy accurately. The only one for which the 3-body recombination rate has been 
calculated is the HFD-B3-FCI1 potential [16]. The scattering length for the HFD-B3-FCI1 
potential is a = 172 oq. This is much larger than either the effective range r s = 14 ao or 
the van der Waals length scale £ v dw = 10.2 ao- The large scattering length a implies that 
4 He atoms have universal properties that are determined by a. These properties can be 
calculated as expansions in powers of the range divided by a. We refer to the leading order 
in this expansion as the scaling limit or the zero-range limit. A thorough analysis of the 
universal properties of 4 He atoms in the scaling limit has been presented by Braaten and 
Hammer [15| . Various 3-body observables for 4 He atoms have been calculated by Platter and 
Phillips to next-to-next-to-leading order in the expansion in powers of the range divided by 
a [13]. The agreement with numerical results from exact solutions of the 3-body Schrodinger 
equation is impressive [H, 3 • 



The 4 He potential supports a diatomic molecule with a single energy level. The binding 
energy of this shallow dimer in the HFD-B3-FCI1 potential is Ed = 1.600 mK [2(|. Experi- 
ence has shown that the universal predictions for 3-body observables are significantly more 
accurate if the dimer binding energy is used as the 2-body input instead of a [15|. Using 
the universal expression for the dimer binding energy in Eq. Q, we obtain the scattering 
length 2 

a He = 164.5 a . (18) 



A convenient conversion constant for 4 He atoms h 2 jm = 43.2788 Koq. 
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The 3-body recombination rate K%(E) for 4 He atoms interacting through the HFD-B3- 
FCI1 potential Jl6j has been calculated as a function of the collision energy from the thresh- 
old to 10 niK [12j. The results were separated into the contributions from J = 0, 1, 2, and 
3. The individual contributions and their sum are shown in Fig. [TJ We can determine the 
value of a*o for 4 He atoms interacting through the HFD-B3-FCI1 potential from the value 
of the 3-body recombination loss rate at threshold: ^(0) = 7.10 x 10~ 28 cm 6 /s. Inserting 
cine = 164.5 a in Eq. (flOl and solving for a„ , we obtain 

a*° = 143.1 a . (19) 

The result is approximately equal to 1.150 au e . The near equality between a^ e and a^e 
reflects the fact that -^3(0) for 4 He is much smaller than the maximum value C majX ha^ e /m = 
3.67 x 10 -26 cm 6 /s allowed by universality, which implies that a/a*o is close to a zero of the 
sine function in the numerator of Eq. ffl2]) . The value of a*o for 4 He can also be determined 
from the binding energy of the excited 4 He trimer, which is E% = 2.62 mK for the HFD- 
B3-FCI1 potential [20j. The resulting value, a*o ~ 146 ao, is consistent with the value in 
Eq. (|T9|) . Its accuracy is limited by the accuracy to which the numerical value of the product 
a*oK* ~ 0.32 is known. 




E[mK] 

FIG. 1: The 3-body recombination rate K%{E) for 4 He atoms (in units of cm 6 /s) as a function 
of the collision energy E (in units of mK) from Ref. [12|]. The four lowest curves at E = 10 mK 
starting from the bottom are the contributions to K%(E) from J = 0, 1, 3 and 2 (thin solid, dashed, 
dotted and dot-dashed curve, respectively). The solid curve at the top is their sum. 



We can use the 4 He results calculated in Ref. 



12| to determine the universal functions 
in the contributions to the 3-body recombination rate for J = 1, 2, and 3. By fitting the 



results of Ref. [12j for J = 1, 2, and 3, the universal functions fj(x) defined by Eq. (fT4"|) can 



be determined in the range < x < 2.5. The behavior of these universal functions near the 
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threshold are 



w 3.43 x 6 , (20a) 
/ 2 (ac) ~ 261 x 4 , (20b) 
/ 3 (x) w 90.3 x 6 . (20c) 

The results of Ref. [12| for J = do not contain enough information to determine all 
the independent universal functions in the expression in Eq. ffToT) that comes from Efimov's 
radial law. However, if we make a few plausible simplifying assumptions, the universal 
functions for J = can be strongly constrained. According to Eq. (j 16 aft , the absolute value 
of the function Su(x) is approximately 0.002 at the threshold x — 0. Our first simplifying 
assumption is that this function remains small compared to 1 and can be neglected for 
x < 2.5. The expression for the 3-body recombination rate in Eq. ffTBI) can then be reduced 
to 

K^(E) ps C max (sin[s ln(a/a* )](l + h 1 (x)) + cos[s ln(a/a*o)}h 2 (x)) 

2 Ha 4 

+i(sm[s \n(a/a*o)}h 3 (x) + cos[s ln(a/a„ )]/i 4 (x)) — , (21) 

where the functions hi(x) are real- valued and vanish at x — 0. In the case of 4 He, the 
coefficients in Eq. fl2TT) are 

sin[s ln(a H e/a"o)] = °- 140 > ( 22a ) 

cos[s ln(a He /a! I o)] = 0.990. (22b) 

A dramatic feature of the J = contribution to K S (E) for 4 He in Fig. [1] is the deep local 
minimum at the collision energy E m - m = 2.1 niK. The rate at the minimum is K^°\E min ) = 
2.3 x 10~ 29 cm 6 /s. The rate decreases by more than an order of magnitude as E goes from 
0.6 niK to 2.1 mK and then it increases by almost an order of magnitude as E goes from 
2.1 mK to 10 mK. This dramatic behavior can be reproduced by slowly varying universal 
functions if the real part of the amplitude inside the absolute value sign in Eq. (121]) changes 
sign at E m i n and if the imaginary part of the amplitude is small compared to the real part 
except near -Emm- Our second simplifying assumption is therefore that the imaginary part 
of that amplitude can be neglected for x < 2.5. Thus the expression for the rate that follows 
from our two simplifying assumptions is 

K^\E) = C max | sin[s ln(a/a* )](l + hi(x)) + cos[soln(a/a* )]/i 2 (£)| 2 /ia 4 /m. (23) 

We expect this to be a good approximation except near those values of a and E at which it 
vanishes, which in the case of 4 He is at E = 2.1 mK. 

The approximation for K^°'(E) in Eq. 0231) st ill involves two unknown functions h\{x) 
and h 2 (x). The J = results for 4 He in Ref. [12J can be used to determine only a specific 
linear combination of those two functions: 

h Hc (x) = 0.140 hi(x) + 0.990 h 2 (x) . (24) 



To determine this linear combination, we fit the J = results of Ref. [12j using the expression 



K${E) = C ma 40.140 + / iHc (x)] 2 K e /m + KS(E min )(2E/( J E; min + E)) 2 . (25) 
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The last term in Eq. (|25l) is a simple model for the term in K^°'(E) that comes from the 
square of the imaginary part of the amplitude in Eq. (T2TJ) . It has the correct value at the 
threshold, where it must vanish, and at E m i n = 2.1 mK. This model affects the fit for hn e (x) 
only in the region very close to E m i n . The resulting function hu e (x) vanishes at x = 0, 
decreases monotonically to —0.135 at x = 1.15, and then decreases further to —0.206 at 
x = 2.5. 



OTHER SYSTEMS WITHOUT DEEP DIMERS 

If the 3-body recombination rate had been calculated for another system of identical 
bosons with a large positive scattering length, such as 4 He atoms with a potential other 
than HFD-B3-FCI1, we would be able to use those results to determine a different linear 
combination of hi(x) and hi(x) than in Eq. (T2"4"|) . We would then know the two independent 
functions hi(x) and h 2 (x), and we would be able to calculate K^°'(E) for any system of 
identical bosons with large positive scattering length. Although we cannot determine hi(x) 
and h,2(x) from the results for 4 He, we can put constraints on them if assume that hn e (x) 
is not unnaturally small because of a near cancellation between the two terms in Eq. (1241) . 
The constraint in Eq. ff24l on h\(x) and h,2{x) can be implemented by expressing them in 
the form 

hi(x) = ^Lh ne (x), (26a) 

h*(z) = ^pM*), (26b) 

where t(x) is an unknown function of x. If t(x) is in the range < t(x) < 1, the two terms 
in Eq. (124"]) contribute to hn e (x) with the same sign. If t(x) is in the range 

-2<t(z)<+3, (27) 

the two terms in Eq. (1251) can have opposite signs, but the absolute value of their sum can 
not be smaller than that of the larger of the two terms by more than a factor of 3. This 
range is consistent with our assumption that hn e {x) is not unnaturally small. Given this 
assumption, we can generate a range of predictions for the 3-body recombination rate for 
other systems with large scattering lengths by varying t(x) in the range given by Eq. (|27[) . 

The 3-body recombination rate K^{E) is shown in Fig. [2] as a function of the collision 
energy E for three values of the ratio a/a*o of the scattering length to the Efimov parameter: 
a/a*o = 1, for which -Kshaiiow(O) = 0, a/a* = 1-150, which corresponds to 4 He, and a/a*o ~ 
4.76, for which i^ s haiiow(0) = C max ha 4 /m. The rate K 3 (E) is the sum of the J = term in 
Eq. (1231) and the J = 1, 2, and 3 terms in Eq. (fl4|) . The rate for each value of a/a*o is shown 
as an error band that corresponds to varying the parameter t defined by Eqs. ( 126]) over the 
range in Eq. (12?]) . Since the error band comes only from the J = term, it always remains 
above the contribution from the J = 2 term. In the 4 He case, a/a*o — 1-150, the band 
reduces to a curve given by the heavy solid line. For a/a* = 1, the rate is dominated by the 
J = 2 term and the error band is always fairly narrow. It extends above the J = 2 curve by 
a factor that is about 1.5 for E < 0.25 Ed and increases to about 3.6 at E — 6.25Ed- For 
a/a*o =4.76, the error band is narrow for E < 10~ 2 E D , but its width increases quickly with 
E. For E > 0.25 Ed, the error band extends above the J = 2 curve by about two orders of 
magnitude. 
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FIG. 2: The 3-body recombination rate K^(E) for atoms with no deep dimers and a large scattering 
length. The rate K%(E) (in units of ha 4 /m) is shown as a function of the collision energy E (in 
units of Ti 2 /(ma 2 )) for various values of the ratio a/a*o'- 1 (solid lines), 1.150 (heavy solid line), 
and 4.76 (dashed lines). The band of predictions for each value of a/a*o corresponds to the range 
-2 < t < +3 for the variable t in Eq. ([26]) . 



Given a range of predictions for K%(E), we can obtain a range of predictions for the 3- 
body recombination event rate constant a(T) given by Eq. Since the results in Ref. 12| 



only allow us to determine the function h^ e (x) in the region x < 2.5, the integral in the 
numerator of Eq. Q can only be evaluated up to the energy E mgx = 6.25 E D . As long as 
T < 0.74 Ed , more than 99% of the weight of the normalizing integral in the denominator of 
Eq. © is in the region x < 2.5. In Fig. [31 a(T) is shown as a function of a for fixed a*o and 
several values of T. The range of predictions is shown by plotting 6 curves corresponding 
to integer values of the variable t defined in Eq. (1261) ranging from —2 to +3. The widths 
of the error bands reduce to zero at the 4 He value a = 1.150 a^o- The error bands remain 
narrow in a region that includes the value a = a* at which K sh£L n ow (E = 0) vanishes. We 
therefore obtain almost unique predictions for a(T) near its local minimum as a function of 
a. The depth of the minimum decreases as the temperature T increases. 



EFFECTS OF DEEP DIMERS 

In this section and in the subsequent section, we consider atoms that have deep dimers. 
In the scaling (or zero-range) limit, the cumulative effect of all the deep dimers on Efimov 
physics can be taken into account rigorously through one additional parameter 77* Q. The 
effects of Efimov physics can be dramatic only if 77* is much less than 1. If the universal 
expression for a scattering amplitude for the case of no deep dimers is known as an analytic 
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FIG. 3: The 3-body recombination event rate constant a(T) for atoms with no deep dimers and a 
large scattering length. The rate constant a(T) (in units of ha^ /m) is shown as a function of a 
for fixed a*o and several values of the temperature T: T = (solid line), 0.007 T*o (dashed line), 
0.07 T*o (dotted line), where ^T»o = ft 2 /(ma 2 ). The 6 curves for each temperature T correspond 
to integer values of t from —2 to +3. 



function of a*o, the corresponding result for a system with deep dimers can be obtained 
without any additional calculation simply by making the substitution 

lna* — ► In a*o — «7*/ s o ■ (28) 

Making this substitution in the factor sin[s ln(a/a*o)] in Eq. ([12]) . the resulting expression 
for the recombination rate at threshold into the shallow dimer is 

K shaacw (E = 0) « C max (sm 2 [s ln(a/a* )] + sinh 2 77*) hcf/m . (29) 

This expression shows that one effect of the deep dimers is to eliminate the zeroes of 
-^shallow (0). Making this substitution in Eq. (1231) . the resulting expression for the J = 
contribution to Shallow is 

K^(E) = C max cosh 2 ^( sin[s ln(a/a* )](l + hi(x)) + cos[s ln(a/a* )]/i2(^)) 

+ sinh 2 77* ( cos[s ln(a/a* )](l + hi(x)) — sin[s ln(a/a Jf0 )]^2(^)) — -(30) 

If there are deep dimers, there is an additional contribution K deep (E) from 3-body re- 
combination rate into deep dimers. A generalization of Efimov's radial law implies that 
K deep (E) must have the form ^ 

K (E) = k\s 31 (x)\\l - e-^) ha 4 

d X 4 |1 - Sll ( X )e 2 ^oln(a/a t0 )-2 J 7*| 2 171 ' 
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where k is the same constant and sn(x) and S3i(x) are the same entries of a symmetric 
3x3 unitary matrix that appear in Eq. (I15p . We can obtain an analytic expression for the 
value at the threshold E = by using the value of sn(0) given in Eq. fll6al) and the limiting 
behavior of s 3 i(x) as x — > 0: 

E = Q) = C max sinh( 27 r g0 )sinh(2^) to* 

4 (sinh (7rs + 77*) + cos 2 [s ln(a/a* )]) w 

Since sinh 2 (7rs + 77*) > 139 is so large, K deep (0) can be approximated with an error of less 
than 1% by omitting the cos 2 term in the denominator of Eq. (|32"1) . To within the same 
accuracy, we can replace sinh(27rso) and sinh(7rso + 77*) by exponentials to get the simpler 
expression 

w^°)~ (1 - e T )Cm " -- (») 

4 m 

One of the simplifying assumption that we used to obtain the expression for K^(E) in 
Eq. (1231) was that |sn(x)| <C 1 for x < 2.5. We can also use this assumption to simplify 
the expression for K deep (E) in Eq. (13~T1) . The only dependence of Kd eep {E) on the collision 
energy E then comes from the factor |si3(x)| 2 in the numerator. The 4 He results do not 
give significant constraints on this factor, because in the expression for K^(E) in Eq. (fl~5l). 
Si3(x) appears in the combination Si2(x)sis(x). We will therefore make the simplifying 
assumption that |si3(x)| 2 is a sufficiently slowly varying function of x for x < 2.5 that we 
can approximate K deep (E) by its value at E — 0, which is given in Eq. fl33l) . 

The dependence of the 3-body recombination event rate constant a(T) on the scattering 
length a is illustrated in Figs. H] and [51 The rate K^{E) in the numerator of Eq. (Q is the 
sum of Shallow (E), for which the J = term is given in Eq. ( |30|) and the J = 1, 2, and 
3 terms are given in Eq. CHI) , and K deep (E), which is approximated by Eq. (I3"3"|) . For each 
value of T and 77*, the range of predictions is shown by plotting 6 curves corresponding to 
integer values of the variable t defined in Eq. ([2"6"j) ranging from —2 to +3. The widths of the 
error bands reduce to zero at the 4 He value a = 1.150 a*o- In Fig. HJ 77* is fixed at the value 
0.01 and a(T) is plotted as a function of a for various values of T. A convenient unit for 
the temperature is T* = h 2 /(kBmal ). As T increases from to 0.007 T^o, the prediction 
for a near a*o is essentially unchanged. As T increases further, the local minimum decreases 
in depth and it tends to disappear around T = 0.1T„o. In Fig. [5l T is fixed at the value 
0.007 T^o and a(T) is plotted as a function of a for various values of 77*. As 77* increases, the 
depth of the local minimum decreases until it disappears around 77* = 0.06. 

Our approximation for the 3-body recombination rate K 3 (E) depends on the fitted func- 
tions hi(x) and h,2(x) in Eq. (I3"U|) and fj(x) in Eq. (TI4T) . If the only energies that contribute 
significantly to the Boltzmann average in Eq. (Q are ones that satisfy E Ed, we can 
use a much simpler approximation. We can approximate K^{E) by the sum of the J = 
term of i^shaiiow(-E') at E = 0, which is given in Eq. ff29l . the leading term in the low-energy 
expansion of the J = 2 term, which is given by Eqs. ([141 and f!20bl) . and i^deep(O), which is 
given in Eq. fl33l) . Inserting these expression into Eq. (Q, we get a simple expression for the 
rate constant for 3-body recombination: 



atotai(T) w 67.1 Uin 2 [s ln(a/a* )] + \ cosh (277,) (1 - e - 2 "*) + 7.8 ( 




(34) 
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FIG. 4: The 3-body recombination event rate constant a(T) for atoms with deep dimers and a 
large scattering length. The rate constant a(T) (in units of ha^ /m) is shown as a function of a 
for 77* = 0.01 and several values of the temperature T: T = (solid line), 0.007 T^o (dashed line), 
0.07 T*o (dotted line), where fceT^o = h 2 /(ma 2 ). The 6 curves for each temperature T correspond 
to integer values of t from —2 to +3. 

APPLICATION TO 133 CS ATOMS 



The Innsbruck group has carried out beautiful measurements of the 3-body recombination 
rate for ultracold 133 Cs atoms 3 in the \f = 3,mj = +3) hyperfine state [8f. By varying the 
magnetic field from to 150 G, they were able to change the scattering length from —2500 a 
through to +1600 a . In this range of magnetic field, the \ f = 3, rrif = +3) state is the 
lowest hyperfine state, so 2-body losses are energetically forbidden. Thus, the dominant 
loss mechanism is 3-body recombination. The van der Waals length scale for Cs atoms is 
(mCe/h 2 ) 1 ^ w 200 a . The range of scattering lengths studied by the Innsbruck group 
includes a region of large negative a and a region of large positive a separated by a region of 
small \a\. In the two regions of large scattering length, few-body physics should be universal. 
An interesting open question is whether there is any relation between the Efimov parameters 

and 77* that characterize the two universal regions. 

In the region of negative a, the Innsbruck group measured the loss rate constant L 3 as 
a function of a at three different temperatures: T = 10 nK, 200 nK, and 250 nK. They 
observed a dramatic enhancement of the loss rate for a near — 850 ao- At T = 10 nK, the 
loss rate as a function of a can be fit rather well by the universal formula for T = in 
Ref. |7j with parameters a'^ = —850(20) a and 77* = 0.06(1). Thus, the large enhancement 
in the loss rate can be explained by the resonant enhancement from an Efimov trimer near 
the 3-atom threshold. More recently, the Innsbruck group has measured the position of 



A convenient conversion constant for 133 Cs atoms is h 2 /m = 1.30339 Kag. 
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FIG. 5: The 3-body recombination event rate constant a(T) for atoms with deep dimers and a 
large scattering length. The rate constant a(T) (in units of ha^ /m) is shown as a function of a 
for temperature 0.007 T*o, where A;_bT*o = h 2 /(mal ), and several values of r/*: (solid line), 0.01 
(dashed line), and 0.06 (dotted line). The 6 curves for each value of 77* correspond to integer values 
of t from —2 to +3. 



the maximum loss rate as a function of the temperature [2l|. Its behavior as a function of 



temperature can be explained at least qualitatively by the dependence of the binding energy 
and width of the Efimov resonance on the scattering length [10, 111 



In the region of positive a, the Innsbruck group measured the loss rate constant L 3 at 
T = 200 nK for values of a ranging up to 1228 a$. Their results in the region < a < 600 ao 
are shown as solid triangles in Fig. [6J The vertical axis is the recombination length p 3 defined 
by 



P3 



( 2m 

W3h 



L, 



1/4 



(35) 



They observed a local minimum in the loss rate for a near 200 clq. This value is near the 
van der Waals length scale (mCg/fr 2 ) 1 / 4 ~ 200 do, so there may be large corrections to the 
universal predictions. The universal prediction for the loss rate at T = is a = K 3 (0)/6, 
where K s (0) is the sum of Eqs. (|29|) and (1331) . By fitting the data for a > 500 ao to this 
expression, they obtained a + = 1060(70) a for the scattering length e 7r ^ 2s °^a^ at which 
the coefficient of a 4 achieves its maximum value. Universality would then imply that the 
minimum should be at a*o = 223(15) ao- The fit was insensitive to the value of 7/* and yielded 
only the upper bound 77* < 0.2. The Innsbruck group also determined the location of the 
minimum directly by measuring the fraction of atoms that were lost after a fixed time. The 
result was a min = 210(10) a , which is consistent with the value obtained by fitting the data 
for a > 500 a . 

We now consider whether our results obtained by scaling results from 4 He atoms can be 
applied to this system of 133 Cs atoms. In the experiments with positive scattering length, 
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the typical peak number density was 5 x 10 13 cm -3 . The corresponding critical temperature 
for Bose-Einstein condensation is T c pa 160 nK. Thus, the temperature T = 200 nK was 
not far above T c . We ignore this complication and calculate thermal averages using the 
Maxwell-Boltzmann distribution as in Eq. (jH]) instead of the Bose-Einstein distribution. 
When a = a min , the universal prediction for the binding energy Ed of the shallow dimer in 
Eq. gives 3 x 10 4 nK, which is about two orders of magnitude higher than the temperature. 
The binding energy E D decreases to about 3 x 10 3 nK at a = 600 a . Thus the temperature 
T = 200 nK is safely in the region T < 0.7 Ed in which the thermal average in Eq. (Q can 
be calculated accurately from scaling the results for 4 He atoms. 

In Fig. ([6]), we compare the universal predictions for the recombination length p% defined 
in Eq. (135]) with the Innsbruck data. We assume ni os t — 3, so that L3 = 3a, and we set 
a*o = 210 do- We plot p3 for T = 200 nK as a function of a for several values of 77*. For each 
value of 7/*, a band of predictions is obtained by plotting 6 curves corresponding to integer 
values of the variable t defined in Eq. ( 12 6 p ranging from —2 to +3. The band has zero width 
at a = 242 a , and the band remains narrow for a near a min = 210 do- The predictions fall 
well below the data in the region a < 300 a , especially if ??* is small enough that there is 
still a local minimum. This suggests that a m \ n is too small for the universal predictions to 
be quantitatively accurate. For a > 500 ao, the scattering length is large enough that we 
expect the universal predictions to be valid. The results of Ref. [§| for a > 400 ao lie within 
the error bands of our universal predictions. Unfortunately, the widths of the error bands 
are too large for a > 500 a to allow us to determine the Efimov parameters a*o and 77*. 
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FIG. 6: The 3-body recombination length for 33 Cs atoms as a function of a for T = 200 nK. 
The data points are from Ref. 0]. The curves are the universal prediction for several values of 
7/*: (solid lines), 0.01 (dashed lines), and 0.06 (dotted lines). The 6 curves for each value of 77* 
correspond to integer values of t from —2 to +3. 
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SUMMARY 



In this work, we have used previously published results on the 3-body recombination 
rate of 4 He atoms as a function of the collision energy to constrain the universal functions 
that govern 3-body recombination processes for atoms with large scattering length. Since 
the scattering length for 4 He atoms differs from the Effmov parameter a*o by only about 
15%, the constraints are very strong near a minimum in the 3-body recombination rate as 
a function of scattering length. We were unable to fit the 133 Cs data from the Innsbruck 
group near the local minimum of the recombination rate by adjusting the Efimov parameters 
a*o and 77*. This is not a surprise, because the local minimum occurs at a value of a that 
is not large compared to the van der Waals length scale. The universal predictions should 
be accurate at larger values of the scattering length for which the recombination rate was 
measured by the Innsbruck group. Unfortunately the constraints from 4 He atoms have 
large error bands at those larger values of a. Our universal predictions obtained by scaling 
results for 4 He atoms can be applied to other systems of identical bosons for which the local 
minimum of the recombination rate occurs at a large value of the scattering length. For such 
a system, our results can be used to give a reliable determination of the Efimov parameters. 

In recent years, an effective field theory has been developed for nonrelativistic few-body 



systems with large scattering length [22|, |23j . It allows the model-independent calculation 
of few-body observables using a systematic expansion in if a. The leading term in this 
expansion is the "universal" term. This effective field theory can be applied to a large 
variety of physical systems ranging from atomic physics to nuclear physics. It has been used 
to calculate the atom-dimer scattering phase shift as a function of the collision energy up to 



the dimer-breakup threshold [15j. The results were used to give universal predictions for the 
resonant dimer relaxation rate at nonzero temperature |24j. The predictions were compared 
with measurements by the Innsbruck group, which revealed a resonant enhancement of the 
inelastic loss rate from a system of ultracold 133 Cs atoms and dimers [251 ] . This effective 
field theory can also be used to calculate the universal functions that determine the 3-body 
recombination rate as a function of the collision energy. Such a calculation would allow the 
Efimov parameters and 77* for 133 Cs atoms to be determined by fitting the Innsbruck data 
for a > 500 ao- It could also be applied to any other system of identical bosons with large 
positive scattering length. Once the Efimov parameters have been determined accurately, 
it will be possible to make quantitative tests of the correlations between different aspects of 
Efimov physics that are predicted by universality. 
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